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I.  INTRODUCTION 


The  objective  of  this  task  was  to  verify  the  exponential  exposure 
model  for  the  vulnerability  of  equipment  to  carbon  fibers  with  contact 
distances  less  than  the  fiber  length  and  to  develop  a mathematical 
model  for  contact  distances  greater  than  one  fiber  length.  The  need 
for  this  work  arises  from  the  fiber  length  spectra  which  have  been 
documented  by  the  TRW  results  from  the  open  air  tests  at  China  Lake. 1 
The  fiber  length  spectra  seen  in  these  tests  have  their  peak  in  the 
region  of  3mm.  If  the  lengths  of  accidentally  released  fibers  are 
short  (<3mm) , it  would  be  beneficial  to  have  a model  which  would  predict 
vulnerability  for  equipment  with  vulnerable  points  greater  than  one 
fiber  length. 


II.  SELECTION  OF  A MODEL 


In  developing  a mathematical  model  to  describe  the  kill  probabilities 
due  to  carbon  fibers,  it  was  found  that  if  a single  fiber  is  capable  of 
effecting  a kill,  that  is  the  gap  between  electrical  contacts  is  less 
than  or  equal  to  one  fiber  length,  then  the  probability  of  kill  due  to 
a single  fiber  is  much  higher  than  the  probability  of  kill  due  to  any 
multiple  fiber  event.  Most  systems  tested  have  been  susceptible  to 
single  fiber  kill  where  the  fibers  were  3.5  to  15mm  in  length. 

The  question  to  be  addressed  here  is;  if  a single  fiber  cannot 
effect  a kill  because  the  gap  width  exceeds  the  fiber  length,  what  model 
should  be  used  to  describe  the  probability  of  kill  due  to  a multiple 
fiber  event? 

Let  F (x)  be  the  distribution  function  of  the  time-to-failure 
random  variable  X,  and  let  f(x)  be  its  probability  density  function. 

Then  the  hazard  rate,  h(x),  is  defined  as 


2.1 


Here  l-F(x)  is  called  the  reliability  at  time  x.  The  hazard  rate,  which 
is  a function  of  time,  has  a probabilistic  interpretation;  namely, 
h(x)  dx  represents  the  probability  that  a device  of  age  x will  fail  in 
the  interval  (x,  x+dx),  or 


h(x)  = lim 


P {a  device  of  age  x will  fail  in  the 
interval  (x,  x+dx) /it  has  survived  up  to  x[ 


Ax-+  0 


Ax 


1 


"Data  Reduction  and  Analysis  of  Graphite  Release  Experiments" 3 Paul 
Liebermanj  et.  al. , TRW  Defense  and  Space  System  Group , 1979. 
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On  the  basis  of  physical  considerations,  h(x)  will  be  chosen  for  a parti- 
cular device.  Given  the  choice  of  h(x)  and  assuming  F(0  ) =0 
and  F (<*>)  = 1,  f(x)  and  F(x)  can  be  determined  as  follows: 


By  definition, 


and 

Now  from  equation  2.1 


f(s)ds  = F(x) 


dF(x) 

dx 


f (x)  . 


dF(x) 

h(x)dx  = l-F(x) 


or 


Thus , 


or 


and 


Taking  the  derivatives 


F(t)  = 1-exp  £-/*  h(x)dx]  2.2 


f(t)  = h(t)  exp  \_-J"  h(x)dx]  . 


The  above  derivation  of  the  distribution  and  density  functions  using  the 
hazard  rate  concept  is  presented  by  Mann,  Schafer,  and  Singpurwalla  in 
Methods  for  Statistical  Analysis  of  Reliability  and  Life  Data. 


The  devices  for  which  we  wish  to  determine  hazard  rates  can  be 
divided  into  categories  based  on  the  number  of  fibers  needed  to  effect 
a kill. 


For  the  case  in  which  a single  fiber  can  effect  a kill,  the  hazard 
rate,  hj(x)dx,  is  the  probability  that  a fiber  arrives  at  a vulnerable 
site  in  the  time  interval  (x,  x+dx) . It  seems  reasonable  to  assume 
that  the  probability  of  a fiber  arriving  in  the  interval  (x,  x+dxl 
would  be  proportional  to  the  concentration  (C)  at  the  site,  indicated  by  the 
position  vector  r,  in  the  time  dx. 
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Let 


hi(x)dx  = a C(r,x)dx 


Substituting  in  equation  2.2,  we  have 


Fi(t)  = 1-exp  [-a  J C(r,x)dx]  . 


0 


f 


By  definition,  the  time  integral  of  concentration  is  exposure  (E) , i.e. 


0 


We  now  have  the  distribution  as  a function  of  exposure,  which 
is  time  dependent. 


Fi(E)  = 1-exp  [-otE(r,x)]  . 


Using  the  method  of  maximum  likelihood, the  best  estimate  of  a can  be  com- 
puted (Appendix  A) . 

For  the  case  in  which  two  fibers  are  needed  to  effect  a kill, 
the  hazard  rate,  h2(x)dx,  is  the  probability  that  the  second  fiber 
arrives  in  the  interval  (x,x+dx)  given  one  fiber  is  at  the  site  at  time 
x.  As  for  a single  fiber  kill,  we  again  assume  that  the  probability  of 
a fiber  arriving  in  the  interval  (x,x+dx)  is  proportional  to  the  con- 
centration in  the  time  interval  dx.  We  are  assuming  that  the  probability 
of  a fiber  being  located  at  the  site  at  time  x is  proportional  to  the 
time  integral  of  concentration,  exposure,  up  to  time  x>  therefore 


h2(x)dx  = BEC  dx 


Substituting  in  equation  2.2,  we  have 


J0 

= 1-exp  [-  BE2] 


It  is  understood  that  concentration  and  exposure  are  both  time  and  position 
dependent.  We  can  again  write  this  distribution  as  a function  of 
exposure,  as 


F2(E)  = 1-exp  [-BE2] 
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and  the  best  estimate  of  3 found  using  the  method  of  maximum  likelihood 
(Appendix  A) . 


This  process  can  be  generalized  to  the  case  requiring  n fibers  to 
effect  a kill.  The  hazard  rate,  h (x)dx,  would  be  the  probability  the 
th  ^ 

n fiber  arrives  in  the  interval  (x,  x+dx)  given  n-1  fibers  located  at 
the  site  at  time  x.  If  the  probability  of  each  of  the  n-1  fibers  located 
at  the  site  at  time  x is  proportional  to  the  exposure  up  to  time  x, 
and  each  fiber  arrives  independently,  then  the  probability  of  n-1  fibers 

being  located  at  the  site  at  time  x is  proportional  to  E(r,t)  n . 

t ll 

The  probability  of  the  n fiber  arriving  in  (x,  x+dx)  is  again  proportional 
to  the  concentration  for  the  period  dx.  Therefore 


h^(x)dx  = y En  1 C dx 


Substituting  in  equation  2.2  and  writing  as  a function  of  exposure, 
we  have 


0 


1-exp  T-YEn] 


The  maximum  likelihood  estimate  can  again  be  found  for  y (Appendix  A) . 

This  distribution  function  is  the  Weibull  distribution  with  shape 
parameter  n and  scale  parameter  y.  For  a single  fiber  kill,  n=l,  this 
gives  a special  case  of  the  Weibull  distribution,  the  exponential  distri- 
bution. 


III.  EXPERIMENTAL  FACILITIES 


A.  Exposure  Chamber 

The  exposure  chamber  used  to  measure  multifiber  vulnerability  was 
a modified  S-280  van.  This  chamber  is  composed  of  two  separate  sections, 
the  actual  test  chamber  and  an  anteroom.  The  fiber  exposure  section 
is  a room  3.5m  x 1.9m  x 2m  with  a window  in  the  side  for  easy  observation 
of  the  experimental  area  during  the  exposure.  The  test  chamber  has 
twelve  muffin  fans  which  are  Variac  controlled,  and  used  to  maintain 
a homogenous  distribution  of  fibers  throughout  the  chamber.  The  use  of 
these  fans  has  produced  a homogenous  concentration  which  varied  between 
103  and  101*  fibers/m3  for  20  minutes  without  the  addition  of  any  fibers. 
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The  low  airflow  allows  the  assumption  that  the  carbon  fiber  velocity  is 
close  to  settling  velocity.  Therefore,  the  experimental  results  are  not 
biased  by  a high  fiber  velocity. 

The  anteroom  (1.9m  high  x 2m  wide  x 2.2m  long)  is  used  for  target 
and  fiber  preparation.  This  room  is  equipped  with  a sticky  paper  pad 
and  a large  vacuum  cleaner  to  prevent  the  migration  of  carbon  fibers 
out  of  the  experimental  area.  For  futher  information  concerning  the 
exposure  chamber,  see  reference  2. 

B*  Fiber  Dispenser 

Fiber  dispensing  is  accomplished  by  an  instrument  designed  and 
fabricated  at  the  BRL3.  This  dispensing  system  is  housed  in  a vertical 
cylindrical  container  20cm  in  diameter  and  lm  high  and  is  located  in 
the  exposure  chamber.  The  exit  port  is  tapered  in  such  a way  as  to 
increase  fiber  velocity  at  that  point  and  eject  single  fibers  into  the 
room.  There  is  a constant  airflow  through  the  dispenser  which  is  a factor 
of  two  greater  than  the  fiber  fall  velocity  (2.5cm/sec).  At  the  base 
of  the  fiber  container  is  an  outlet  which  emits  a short  burst  of  high 
pressure  air  at  regular  increments.  Both  the  air  burst  duration  and 
repetition  are  electronically  controlled.  The  short  burst  of  air  serves 
to  lift  a large  clump  of  fibers.  The  single  fibers,  because  of  the 
steady  upward  flow  of  air,  will  continue  up  the  container  and  be 
ejected.  The  clumps,  because  of  their  greater  fall  velocity,  will 
settle  back  to  the  bottom  and  be  reelevated  by  one  of  the  following 
bursts.  The  fibers  used  in  the  dispenser  are  precut.  The  fibers  dis- 
pensed are  90  per  cent  single  fibers  with  no  noticeable  length  breakup. 
Approximately  lgm  of  precut  fiber  will  produce  an  exposure  of  1 to  2 x 
10'  fiber-sec/cubic  meter  (f*  s/m3)  of  7.5mm  long  fibers  in  the  exposure 
chamber  during  a 30  minute  trial . 

C.  Fiber  Detection 

There  were  two  methods  of  fiber  detection  used  during  this  task. 

The  active  method,  which  was  the  BRL  ball  gauge,  was  used  to  measure 
the  exposure  at  the  time  of  failure.  The  ball  gauge  uses  electronic 
circuitry  to  count  the  number  of  fibers  present  which  are  greater  than 
one  half  the  nominal  length  being  dispensed.  The  data  was  stored  in 
a multichannel  analyzer  in  the  multichannel  scaling  mode.  Thus,  the 
record  was  the  concentration  as  a function  of  time,  the  sum  of  which  is 
exposure. 

The  passive  method  of  detection  was  a sticky  paper  sample  39mm 
square  which  was  used  to  measure  the  fiber  length  sprectrum  dispensed. 
Because  the  purpose  of  the  work  was  to  express  the  probability  of 
failure  as  a function  of  exposure  with  varying  fiber  lengths, 


A.  Croce  and  G.  A.  Dunn , " Filter  Evaluations  - Project  HAVE  NAME",  ECOM 
Research  and  Development  Technical  Report  No.  ECOM- 42 86,  January  1975 , 
SECRET. 

Private  Communication,  Neil  Wolfe  of  BRL. 
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the  concentration  and  the  fiber  length  had  to  be  carefully  monitored.  A 
typical  concentration  profile  can  be  seen  in  Figure  1 and  a fiber  length 
spectrum  in  Figure  2.  For  further  information  about  the  detectors,  see 
reference  4. 


IV.  TEST  DESCRIPTION 


A.  Test  Objective 

The  objective  of  a test  was  to  measure  simultaneously  the  exposure 
at  the  time  of  the  failure  for  the  different  electrode  gaps.  These 
tests  were  performed  using  7.5mm  Hercules  HMS  fibers.  The  fiber  vel- 
ocity was  settling  velocity  (2.5cm/sec),  and  the  targets  were  oriented 
at  0,  45,  and  90  degrees  to  the  horizontal. 

B.  Target  Description 

1,  Five  Space  Board 

The  contact  spacing  of  interest  were  L/3,  L/2,  L,  2L,  and  4L,  where 
L is  the  nominal  fiber  length  being  dispensed  (7.5mm).  The  five  different 
space  targets  were  located  on  a phenolic  board,  16cm  wide  and  25cm  long. 
Figure  3 shows  the  board.  Each  target  consists  of  two  strands  of  20 
gauge  wire,  5cm  long  separated  by  the  required  distance.  To  insure 
no  interaction  between  targets,  a wooden  insulator,  1.25cm  high,  was 
used  to  separate  targets.  During  each  trial  there  were  four  boards  placed 
in  the  exposure  chamber. 

2.  Three  Space  Board 

During  the  first  series  of  tests,  the  data  showed  an  exposure  limit 
for  the  large  electrode  gaps.  To  investigate  this  phenomenon,  it  was 
decided  to  construct  a target  board  with  very  large  electrode  gap  to 
fiber  length  ratios,  5L,  10L,  and  20L.  Figure  4 is  a photograph  of  this 
board.  There  are  no  insulating  barriers  between  targets  on  these  boards 
because  it  was  constructed  in  such  a way  that  the  adjoining  wire  elements 
were  at  the  same  electrical  potential,  and  thus,  there  could  be  no 
possibility  of  any  target  interaction.  These  targets  were  constructed  the 
same,  two  5cm  lengths  of  20  gauge  wire.  Because  these  tests  were  only 
to  check  the  previous  data,  the  trials  were  only  performed  at  one  angle, 
0°. 

C.  Data  Acquisition 

As  was  mentioned  earlier,  an  individual  target  consisted  of  two 
5cm  long  strands  of  20  gauge  wire  with  the  appropriate  electrode  gap. 

One  electrode  was  connected  to  a lOOVdc  source  and  the  other  to  ground 
through  a 9.1ft  resistor.  When  the  electrode  gap  was  bridged,  whether 
by  one  or  many  fibers,  a strip  chart  recorder  would  measure  the  voltage 
drop  across  the  resistor.  Figure  5 is  a block  diagram  of  the  measurement 
circuit . 


4 

John  A.  Morrissey > W.  I.  Brannany  S.  C.  Thompson , "Calibration  of  BRL 
Ball  and  Sticky  Cylinder  Detector  Systems" , Ballistic  Research  Laboratory 
Technical  Report , ARBRL-TR-02079 , June  1978.  (AD  #B029204L) 
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Figure  1.  Typical  Concentration  and  Exposure  Plot 
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Figure  2.  Typical  Fiber  Length  Spectrum 
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Figure  4.  Photograph  of  Three  Target  Board 


+ — 

100  Vdc 

INPUT 


RECORDER 

* 

f f 

1 j AMP 

rvo— 


9.1  a 

A/W- 


9.1  a 

AAA/ — 


TARGET 

ELEMENT 


19 


Figure  5.  Target  Monitoring  Circuit 


D.  Test  Procedure 


After  the  target  boards  had  been  placed  in  the  chamber,  the  carbon 
fiber  dispensing  began.  The  fibers  used  were  7.5mm  Hercules  HMS, 
which  has  a conductivity  of  3000  ohms/cm.  A data  event  was  defined  as 
a noticeable  current  spike  ( _>  2ma)  on  both  the  voltage  and  ground 
sides  of  a target.  The  exposure  was  continued  until  there  were  two 
data  events  at  each  target,  at  which  time  the  exposure  was  terminated. 

V.  DATA  ANALYSIS 


The  exposures  to  failure  for  each  test  are  given  in  Tables  I thru 
VI.  The  empirical  distribution  and  theoretical  distribution,  as  derived 
in  Section  II,  are  plotted  on  Figures  6 thru  23.  The  single  fiber 
model  is  used  for  L/3,  L/2,  and  L.  For  2L,  4L,  5L,  and  10L  the  minimum 
number  of  fibers  needed  to  bridge  the  gap  is  the  number  of  fibers  used 
in  the  model.  A fourteen-fiber  model  is  used  for  20L  since  anything 
larger  overflows  the  computer.  The  figures  give  1)  the  orientation  of 
the  target  - the  contact  spacing  in  fiber  lengths,  2)  the  maximum 
likelihood  estimate  of  the  parameter  y,  (or  a or  8)  3)  the  number 

of  tests  and,  4)  the  Kolmogorov-Smirnov  goodness-of-fit  statistic  T.  5 
An  asterisk  beside  the  T statistic  indicates  significance  at  the  a=.05 
confidence  level. 

Thirteen  of  the  eighteen  hypothesized  distributions  fit  the  data 
using  the  Kolmogorov-Smirnov  goodness-of-fit  test  at  the  95%  confidence 
level.  The  five  distributions  which  do  not  fit  are  from  the  data 
sets6L/3,  L/2,  and  4L  at  a 45°  orientation,  10L  and  20L.  The  model  that 
best  fits  the  data  45°-L/3  and  45°-L/2  is  a model  using  half  a fiber 
(n=.5),  shown  on  Figures  24  and  25.  It  is  possible  that  the  fiber  being 
two  or  three  times  longer  than  the  gap  enhances  the  kill  probability. 
However,  the  four  other  tests  at  the  same  fiber  lengths  do  not  indicate 
an  enhanced  kill  probability. 

The  data  set  45°-4L  is  best  fit  using  a single  fiber  model.  Figure  26. 
Since  a 7.5mm  fiber  length  and  a 30mm  gap  make  this  impossible,  no 
explanation  is  offered.  Figure  27  shows  the  family  of  curves  for  the 
number  of  fibers,  n = 1,2, 3, 4, 5 and  10. 

The  last  three  sets  of  tests  having  gap  widths  5L,  10L,  and  20L, 
were  conducted  to  see  if  there  was  an  exposure  limit  for  the  large 
electrode  gaps.  The  model  that  fits  each  of  these  data  sets  best 
is  the  five  fiber  model.  Figures  21,  28,  and  29.  This  would  indicate 
that  a gap  width  greater  than  or  equal  to  five  times  the  fiber  length 
needs  an  exposure  that  essentially  covers  the  target.  The  five  fiber 
model  would  then  be  the  limiting  distribution  to  describe  the  situations 

** Appendix  B gives  the  distribution  of  the  K-S  statistic  for  the  Weibull 
distribution  with  scale  parameter  unknown. 

° The  best  fit  for  a set  of  data  is  the  model  having  the  smallest  T value. 
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TABLE  I.  Summary  of  L/3  Electrode  Gap  Data 


EXPOSURE  (f-s/m3) 

TARGET  ORIENTATION 

0°  45°  90° 


Test  No. 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


1.19xl05 
3.69xl05 
3.67x10’’ 
1.68xl05 
5.87xl05 
5 . 95xl04 
7.17xl05 
1.14xl06 
4. 09x1 05 
6. 7xl04 
9. 9xl04 
4.0xl03 


1.03xl05 
1 . 53xl05 
7 . 95xl04 
2.0xl03 
2.8xl06 
1.35X105 
6 . 5xl03 
6. 5xl03 


1.32xl06 

3.29xlOS 

1.48xl06 

1.62xl06 

1.02xl06 

4.79xlOS 

8.05xl0S 

5.52xlOS 


TABLE  II.  Summary  of  L/3  Electrode  Gap  Data 


EXPOSURE  (f-s/m3) 


TARGET  ORIENTATION 
0°  45°  90° 

Test  No. 


1 

4 . 35xl05 

5.40xl06 

1.12x10 

2 

2.54xl05 

1 . 02xl06 

2.19x10 

3 

1.82xl06 

6.5xl03 

1.09x10 

4 

1.68xl05 

9.4xl04 

1.02x10 

5 

5 . 87xl05 

1.13xlOS 

2.77x10 

6 

5 . 95xl04 

5. 6x1 04 

1.71x10' 

7 

2 . 15xl04 

3. 5xl02 

1.47x10' 

8 

4.09xl05 

6. 85xl04 

9 

4.52xlOS 

10 

1.99xlOS 

11 

2.6xlOS 

12 

1 . 94xl06 
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TABLE  III.  Summary  of  L Electrode  Gap  Data 


EXPOSURE  (f-s/m3) 

TARGET  ORIENTATION 
0°  45°  90° 


Test  No. 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


4 . 18xl06 
4 . 28xl05 
9.58xl05 
9.79xl05 
1 . 52xl06 
7.76xl05 
2.54xlOS 
1.45xl06 
8.08xl0S 
6.02xl06 
2 . 98xl06 
1.26xl06 


6.75xlOS 
1.13xlOS 
1.44xl06 
8.27xl0S 
6.05xl06 
4 . 10xl06 
1.63xl06 
9.4xl04 


1.20X107 
1.82xl06 
1.43X107 
6.39X106 
5.33X106 
4. 70x1 06 
8.87xl06 
7.42xl06 
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TABLE  IV.  Summary  of  2L  Electrode  Gap  Data 


EXPOSURE  (f-s/m3) 

TARGET 

ORIENTATION 

Test  No. 

0° 

45° 

\o 

o 

o 

1 

1.62xl06 

9.01X106 

2 . lxlO7 

2 

8.19xl06 

4.34xl06 

4.03x10 

3 

1 .71xl06 

2.85xl06 

3.03x10 

4 

4 . llxlO6 

6.47xl06 

1.55xl07 

5 

5.21xl06 

7.84xl06 

2.39xl07 

6 

3.48xl06 

4.89xl06 

3.76xl07 

7 

4 . 74xl06 

3.61xl06 

8. 9xl07 

8 

5.25xl06 

5. 46x1 06 

7.70xl07 

9 

7.80106 

10 

3.38xl06 

11 

1 .82xl06 

12 

2. 66x1 06 
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TABLE  V.  Summary  of  4L  Electrode  Gap  Data 


EXPOSURE  (f-s/m3) 


TARGET 

ORIENTATION 

0° 

45° 

90° 

Test 

No. 

1 

6. 03x1 06 

1.67xl07 

4. 03x10 

2 

1.74xl06 

1 . 09xl07 

5.05x10 

3 

4.83xl06 

2.85xl06 

2.39x10 

4 

9. 39x1 06 

6.87xl06 

3.76x10 

5 

4.38xl06 

2.22xl07 

8.90x10 

6 

5.38xl06 

2.15xl07 

7.70x10 

7 

6.11xl06 

6.23xl06 

8 

6.33xl06 

5.30xl06 

9 

6. llxlO6 

10 

6. 40x1 0^ 

11 

4. 84x1 06 

12 

3. 38xl06 
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TABLE  VI.  Summary  of  Large  Electrode  Gap  Data 


EXPOSURE  (f-s/m3) 

ELECTRODE  GAP  SPACINGS 
5L  10L  20L 


Test  No. 


1 

7.18xl06 

6.25xl06 

6.25x10 

2 

7.59xl06 

6.89xl06 

6.03x10 

3 

9.63xl06 

6.25xl06 

5.89x10 

4 

9.44xl06 

6.63xl06 

6.81x10 

5 

5.46xl06 

5.65xl06 

5.65x10 

6 

8 . 62xl06 

7. I2xl06 

5.77x10 

7 

1 .10xl07 

8.97xl06 

8.12x10 

8 

1 .02xl07 

8. 97xl06 

7. 36x10 

9 

9 . 07xl06 

6.18xl06 

5.28x10 

10 

8.87xl06 

8.77xl06 

7.34x10 

11 

9.36xl06 

5.78xl06 

5. 78x10 

12 

8.97xl06 

9.55xl06 

9.24x10 

13 

7. 92xl06 

9. 36x1 06 

6.18x10 

14 

8.77xl06 

9.34xl06 

l.OxlO7 

6 

i 

i6 

6 

6 

,6 

6 

I 

,6 

,6 

,6 

6 

6 

6 

,6 
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Figure  6.  Plot  of  the  O-L/3  Data  and  the  Theoretical  Distribution 

For  a 1 Fiber  Model 
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Figure.  7.  Plot  of  the  45-L/3  Data  and  the  Theoretical  Distribution 

For  a 1 Fiber  Model 
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Figure  8.  Plot  of  90-L/3  Data  and  the  Theoretical  Distribution 

For  a 1 Fiber  Model 
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Figure  9.  Plot  of  O-L/2  Data  and  the  Theoretical  Distribution 

For  a 1 Fiber  Model 
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Figure  10.  Plot  of  45-L/2  Data  and  the  Theoretical  Distribution 

For  a 1 Fiber  Model 
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Figure  11.  Plot  of  90-L/2  Data  and  the  Theoretical  Distribution 

For  a 1 Fiber  Model 
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Figure  12.  Plot  of  0-L  Data  and  the  Theoretical  Distribution 

For  a 1 Fiber  Model 
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Figure  13.  Plot  of  45-L  Data  and  the  Theoretical  Distribution 

For  a 1 Fiber  Model 
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Figure  14.  Plot  of  90- L Data  and  the  Theoretical  Distribution 

For  a 1 Fiber  Model 
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Figure  15.  Plot  of  0-2L  Data  and  the  Theoretical  Distribution 

For  a 2 Fiber  Model 


Figure  16.  Plot  of  45-2L  Data  and  the  Theoretical  Distribution 

For  a 2 Fiber  Model 
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Figure  17.  Plot  of  90-2L  Data  and  the  Theoretical  Distribution 

For  a 2 Fiber  Model 
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Figure  18.  Plot  of  0-4L  Data  and  the  Theoretical  Distribution 

For  a 4 Fiber  Model 
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Figure  19.  Plot  of  45-4L  Data  and  the  Theoretical  Distribution 

For  a 4 Fiber  Model 
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Figure  20.  Plot  of  90-4L  Data  and  the  Theoretical  Distribution 

For  a 4 Fiber  Model 
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Figure  21.  Plot  of  0-5L  Data  and  the  Theoretical  Distribution 

For  a 5 Fiber  Model 
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Figure  22.  Plot  of  0-10L  Data  and  the  Theoretical  Distribution 

For  a 10  Fiber  Model 
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Figure  23.  Plot  of  0-20L  Data  and  the  Theoretical  Distribution 

For  a 14  Fiber  Model 
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Figure  24.  Plot  of  the  45-L/3  Data  and  the  Theoretical  Distribution 

For  a 0.5  Fiber  Model 
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Figure  25.  Plot  of  the  45-L/2  Data  and  the  Theoretical  Distribution 

For  a 0.5  Fiber  Model 
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Figure  26.  Plot  of  the  45-4L  Data  and  the  Theoretical  Distribution 

For  a 1 Fiber  Model 
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Figure  27.  Plot  of  the  45-4L  Data  and  the  Family  of  Theoretical  Distributions 
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Figure  28.  Plot  of  the  0-10L  Data  and  the  Theoretical  Distribution 

For  a 5 Fiber  Model 
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Figure  29.  Plot  of  the  0-20L  Data  and  the  Theoretical  Distribution 

For  a 5 Fiber  Model 


having  a large  (i.e.  _>  5)  gap  width  to  fiber  length  ratio. 


VI.  CONCLUSIONS 


The  Weibull  distribution 

F (E)  = l-e"yE 

where  n is  the  number  of  fibers  needed  to  bridge  the  electrode  gap,  and 
y is  estimated  from  the  data  using  the  method  of  maximum  likelihood  is 
a very  acceptable  and  even  desirable  model  because  of  its  simplicity 
in  describing  failures  due  to  carbon  fibers.  When  n = 1,  we  have  an 
exponential  distribution  which  is  widely  accepted  in  the  case  of 
a single  fiber  kill. 

From  the  tests  with  gap  width  greater  than  or  equal  to  five  times 
the  fiber  length,  the  limiting  distribution  appears  to  be  a five  fiber 
model.  That  is,  the  situation  in  which  a target  is  essentially  covered 
by  fibers  in  order  to  effect  a kill  is  best  fit  by  a five  fiber  model. 
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APPENDIX  A 

Maximum  Likelihood  Estimate  of  the  Scale  Parameter  for  a Two 
Parameter  Weibull  Distribution 


The  likelihood  function,  L,  is  the  likelihood  that  a random  variable 
E assumes  particular  values  Ej,  E2,  ...  , Er,  that  is,  the  probability 

of  obtaining  a sample  with  outcomes  (exposures  to  failure)  Ej , E2,  .... 

Er»  In  particular,  if  Ej,  E2,  ...»  Er  is  a random  sample  from  the 

density  of  f(E;y),  then  the  likelihood  function  is 

L = f (Ei ; y)  f (E2;  Y)  ...  f (Er;  y) . 

In  our  tests  we  have  two  possible  outcomes;  1)  a failure  is  observed 
at  exposure  E.  or  2)  no  failure  occurs  and  the  test  is  terminated  at 
exposure  E^. 

The  probability  to  failure  distribution  in  which  we  are  interested 
is 


Fn(E.)  = l-exp(-YEin) 


The  density  function  associated  with  a failure  is  therefore 

fn  = nYEin  l exP  C-yEin) . 


The  density  function  associated  with  a test  terminated  (survival)  at 

E.  is 
1 


fn  (E.)  = exp  (-YE.n). 

Therefore,  if  in  a sample  of  r tests  there  are  k failures,  the  likelihood 
function  is 


L 


(E.n"^exp(-Y  E E.n) 
(E.n"^exp(-Y  ^ E.n). 


exp 


(-Y  2 
i=k+l 


3 
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Maximizing  the  natural  log  of  the  likelihood  function  is  equivalent 
to  maximizing  the  likelihood  function  itself,  so  that 

^ k 1 r 

<*£  = In  L = k (In  n)  +k  (In  y)  + E (In  E.n  - y E (E.n) 


i=l 


i=l 


= £ . I fE 

dY  Y V 1 } 


= * - I (E.n) 

Y i=l 

= k/  Z (E.n) 
Xi  = l 1 


Y is  the  maximum  likelihood  estimate  of  y and 
mean  exposure  to  failure. 


^ / y is  called  the 
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APPENDIX  B 


Distribution  of  the  Kolmolgorov-Smirnov  Statistic  for  the  Weibull 

Distribution  with  Scale  Parameter  Unknown 


The  Weibull  distribution  derived  in  Section  II  is 

F(E)  = l-exp[-yEn] . 

For  a given  n,  and  y estimated  from  the  data  using  the  method  of  maximum 
likelihood,  the  distribution  of  the  Kolmolgorov-Smirnov  (K-S)  Statistic 
will  be  determined. 

A set  of  uniform  random  numbers,  R.  , are  generated  and  solving 
the  following  equation  a set  of  exposures  to  failure,  E.  , are  det- 
ermined for  a given  n and  y 

= exp[-yEin] 

/ 1 y/n 

and  Ei  = ( - In  ^ ) B.l 


Since  y is  unknown,  it  is  necessary  to  estimate  y from  the  data. 

Using  the  method  of  maximum  likelihood  (Appendix  A) 

Y = k/  E E n 
i=l 

where  r is  the  number  of  tests  and  k is  the  number  of  failures.  Sub- 
stituting E^  from  Eq.  B.l,  we  have 

r 

y = -yk/ln  Ri 

Therefore, 

F(E)  = 1-exp  [-yEnJ 

r 

= l-exp  [ykEn/ln  11  Ri] 


The  K-S  statistic,  T,  is  defined  as  follows: 


T. 

1 


Is  (E(i))  - F (E(i))| 
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Ti-  ■ Is  <E(i-i)>  - p <EtiP 


and  T = max  jl\,  Ti'| 
i=l, . . . ,k 

where  S(E^)  is  the  empirical  distribution  function  and  S(Eq)  = 0 
by  definition.  The  subscript  (i)  indicates  the  ith  order  statistic. 

Note:  R. <1 , E.  > 0,  and  In  R.<  0.  E < E0  < . . . < E,  and 
i i i 12  k 

E^  <Ej  implies  R^  > R^  so  that  the  rankings  on  the  E^fs  and 

R.fs  are  in  reverse  order. 
i 

Ti  ■ I s <E(iP  - F (E(i))  I 


exp  [ykCE^p  /In 


Similarly, 


n- 

— ^ — -+  exp  JykCE^p  / 


In 


T.  - = 
1 


exp  [vU.B'Ji,)  /In  ( ^ Rj) 


P|T>^T*}  = P |max  { Ti,  Tt-  | > T*| 


= 1 - P)T.<  T*  H T.<  T*.  V i.i- 


p { V T*}  - p { |*r  * ^kEn(i) /ln  (."  |iT*} 

Solving  for  E^,  we  have  P i T\  £ T*|  = 

r (in  n Ri)[m  (¥  -T.31,/n>  «...  ki-nUf.J 

L i=l  yk  J UJ  L i=l  yk  J 

Substituting  for  E ^ from  Eq.  B.l  and  using  the  fact  that  the  E^s 

and  R.fs  are  in  reverse  order. 
i 


1/n 
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jTi  iT*j  * 


r R.  k'i 


(In  n j)  [In  ( k 
' i=l  yk 


T*)] 


1/n 


(-i  111  R(k*l-i)) 


1/n 


^ R ^ 

(In  n 1)  [In  ( k + T*)  ] 
i=l 


Therefore, 

P 


> 

{ Ti  1 T*}  = p j IllL  ,n  Ri)  [ln( 


yk 


k-i 


1/n  J 


•k  - T*)  ] > -In  R..  . ..  > 
■ 1 - (k+l-i)  - 


Using  a similar  argument 


(in  h *i,  [ln(¥.  T.)]i 
i=l  k ' 

{Ti_lT*}  = P | _(ln  n Rj)  [in  C k"  - T*)]  > 

(In  S Ri)  [lnfnH-  + T*)]} 
i=l  k 1 


-In  R,.  , . > 
(k+l-i)  — 


Therefore,  the  distribution  of  the  T statistic  of  the  Weibull  dist- 
ribution with  the  scale  parameter  estimated  from  the  data  is  independent 
of  the  shape  parameter,  n,  and  the  tables  of  the  distribution  of  the  T 
statistic  for  the  exponential  distribution  with  mean  unknown  6 are 
appropriate. 


6 


Hubert  W.  Lillie fors,  "On  the  Kolmogorov- Smirnov  Test  for  the 
Exponential  Distribution  with  Mean  Unknown",  JASA , March, 1969,  pp.  387-389. 
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